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ABSTRACT 


This thesis gives a number of approximations and bounds 
for the renewal function im an ordinary renewal process. 
Each approximation and bound is calculated for the uniforn, 
gamma and hyperexponential distributions and compared with 
the renewal function for these cases. They are also calcu- 
imated for the log-normal distribution an@ compared) wi theue. 
sults of the simulation of the renewal function. Results 


are tabulated and shown graphically. 
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I. INTRODUCTION 


Let {X, ,X5,...} be a sequence of independent identi- 
cally distributed non-negative random variables. Let the 


mestribution funetion of Xs, 1 =6l, 2s FOS ret 


the nth partial sum, and So = 0. Define the random vari- 
able 


(Gey) = Sroyn a 2) Se Se I 


n 
which iS the number of renewals occurring before and in- 


eluding t. The renewal function, M(t), is defined as 
Mie) = EGG). ties (1) 


Let F(t) be the n-fotd convolution: of F with 1¢sele eae 
melationship between M(t) and F(t) is given by the follow- 


ing, 


M(t) =) FA(t), 
ined | 


which leads to the well-known renewal equation, 


t 
M(t) = F(t) “7 Meee (ca, (2) 
O 
The purpose of this thesis is to investigate approxima- 
tions and bounds for the renewal function, since Equations 
(ie and (2) can only be Solved for Special eases Ob pune arc. 


Errol lone Tune: Tom —. 








In Section I] formulae are collected for various ap- 
proximations and bolnds Senat Nave sappearcdalim tie litera. 
ture (Barlow §& Proschan (1965), Bartholomew (1963), and 
Butterworth and Marshall (1971)). These formulae are 
applied to the uniform, gamma, hyperexponential and log- 
moral distri bubronseIneoccelOne li tie —mesul ts 10d scom- 
putations aréespresented in both tabular and graphical 


morms fOr each approximation applied to each distribution. 








If. APPROXIMATION FORMULAE 


Throughout this thesis the following notation will be 








used: 
1 
E[X; ] ~ Oe ) 
ME[X=] 
k = 5 - 1, 

and - 

F(t) - af (Cee aemeerGU) ) CU aetna met 

O 


With this notation De 1S an Caui librium excess ars tuition 
weed renewal process, and from Smith (1958), 1f F as non- 


teat Lice , 


are “ pak Ny AS e Us 8 Go) ON 
Mit j Ne 7 tee ee) e eo. 


ibis Equation, together with (2), plays an important pant 
mieche development of approximations to M(t): 
1. <A Simple Lower Bound 
Butterworth and Marshall (1971) show that for any 
renewal process 
n n 
WiC ej] ee wee es ». F(t) - yy Eee Fi4 Ct) - F(t), 
k=] k=1 
ees On (4) 
miere F(t) = 1, and when n=0 both summations are assumed” to 
pe empty. 8lhe right-hand sidewot (4) gives an increasing 
sequence of lower bounds which converge monotonically on 


M(t). The case n=l will be used in this thesis and we define 


Ge =e At F(t), 6D fil. (5) 





= 








2. An Improved Lower Bound 


The following general lower bound was found by 


Barlow and Proschan (1965). Define 


At 
Fey i (6) 


B(t) = 





Minen B(t) < M(t). 9 Sance Bit) = (A(t) /F,(t)), 1t isMcileanr 
that B is an improvement over A when F(t) ele 
De pve wel sunaSse ite siats  NiRUE 
A distribution of a non-negative random variable is 


called NBUE (new better than used in expectation) if 


mor all t where F(t} < 1. Equivaientiy F is NBURRai fr er see, 
ment) for all t. For détails concerning this and other 
classes of distributions frequently used in reliability 
theory see Marshall and Proschan (1970). Butterworth and 


Pesshall (1971) show that 1£ F 1s NBUE, 


n n 
M(t) < At + ) A(t) - ) FoF y(t) , m2 0, (8) 
k=] k=1 
where again Ges = 1, and when n=0, both summations are 


taken to be empty. The right-hand-side of (6) gives a de- 
creasing sequence of upper bounds which converge monotonic- 


ally to M(t).! The case n=2 will be used in this thesis. 


1 

If the mequalitty 1s reversed in“(/), Fears called Niue 
(new worss than used in expectaticn). In this case (8) gives 
an improved sequence of lower bounds over the sequence given 


Dye 4 le 





We define 


C(t) = At + F(t) - FL(t) (9) 


Atel F(t)e 
4. Upper Bound When F is IFR 
A distribution is said to have increasing failure 
mace {15 IER) af Log(l-BGt)) rms concave in t..) Under thiceas- 


sumption Barlow §& Proschan (1965) showed that if we define 


. AG) . (10) 


D(t) = 


men D(t) > M{t). 
5. Approximations for M(t) 
Bartholomew (1963) derived the following approxima- 
tion for the renewal density m(t} = (aM{cj}/edt} whem F has 


density f; 


Ge) = 


m(t) iS 16 (76) a Ea) 


The following integral of this approximation is used as an 


approximation for M(t), 


t 
E(t) = F(t) +A J ak du. (11) 
O 


Another form of the renewal function can be derived from (1) 


or (2). One can Showsethat 


ic 
M(t) = At *f [1 - F,(t-u)] dM(u) - F(t). (12) 
i | 


If we approximate dM(u) on the right hand side by Adu we get 


MOLLere approximation £or M. Thus define, 


10 





re 
G(t) = at + a J [1 - F,(u)] du - F(t). (13) 


O 

Recall that A(t) and B(t) gave simple lower and upper bounds 
on M for F restricted to NBUE. However, (A(t)-At)7>-1 and 
(B(t)-At)-0 as t20. From (3) for F non-lattice Ve 1s known 


that (M(t)-At)~k, so define 


H(t) = At - F,(t) t (lle): eer) (14) 


the convex combination of A and B with the correct asymptotic 
behavior. 
6. Distribution Tested 
The following distributions with specified param- 


eters are used to compare the bounds and approximations for 


M(t). 
a. Uniform Distribution (A=1) 
o . 
CE) eee a Oia << 72 
= ] Lt Zee 
Then 
te. 
Ge el a oa We OP ee Ss 
= ] ie te 
These give: 
te | : 
A(t) = “= Tae 10) 8g ge See 
= 67) ok sie Cs 
Go ee ie ee Se 
whoie - 
=t- 1 iste Se 


i 








CGE) 
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E(t) 


G(t) 


Bate) 


For this case, 


M(t) 


t= t 

A 25 a sae 

1e 1 
2t : 

q-t ee 

t if 
4 2 : 

4log ( q-t ) zi > lisa 
1 : 

tp 446g) > re if 
te Ss 

TE = 7 2 | © ae 
mat 

t z te 
tet 

om v = The 
1 : 

cae iz ite 


jel tu 
yap (n-M ec” [-1 if 
n=0 


b. Gamma Distribution (A=1) 


nile) 


F,(t) 
These give: 


A(t) 


B(t) 


it 


ies Gl 22) e SF 


cer lett) 5 BU 
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AG ae ee a. 


c ‘ 
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——————_- -l, 
1 - (1 + t) e 


(j=1) (2 = ee 
j = 1,2,... 





C(t) 


D(t) 


E(t) 


Gree) 
H(t) 
For this case 


M(t) 


Cc. 


F(t) 


F gift) 


These give: 


A(t) 


Bee 


C(t) 


D(t) 


ll 


Wee ge a 
l1-e one 2S i 
at 
-2u - 2 
ee ee f Gl xe eC a 
A 1 -e : - 
1 1 -2t 1 
Co ee - 7 
1 1 -2t 1 
t 7 (t x7) e “as 
1 -4t 1 
oe © q° 


Hyperexponential Distribution (A=1) 


Zi 
3 _ =2t 1 rece 
7 (1 € era € Ne 
2, 
5 26 5 5 
g (1 c ee an c ve 
Z 
-— t 
cas 3 ou + = e : =a, 
St a 
ze 
3(1 - 2 & ) + S(l -e : ) 
3 
3 aaa g 5 
Bo el 1 ag GL ome 
an 
Zepocl. eae eo kee : yl 
-E t 
SC 1 ace egese 





2 
= ie 
Ce tees cy See GY} 
q q 
2 
t -=— U 
oY 
+f ace) el a 
O = 20) ¢ 3 
OME 5) FanGr 6c ) 
_— 
3 -2t, . 15 5 
Glee) S bes tm ilo = «@ as ne ik =e ie 
2, 
15 -2t 3 
H(t) = t+ 72(1-e 4) -a(1l-e? ). 
For this case 
4 
3 can 
C2) ees 7p Cee ) 





Let 
| t _X* 
o(t) = : f e ’ dx 
V27 
Then 
F(t) = 6 (log t) 
i 
F(t) = te * (1 - @(log t)) + ® ((log t) - 1). 


maese give: 


i 
A(t) =e é Ge willoe 2) = OO LL@e wi) 2 Il. 
=a 
u 
B(t) = a i ie 


57 Beno Helos 2) 26 Clos 4) 08 





i 


C(t) = (6 © 4 Wiiclos teen (log tine 


1 
a 
D(t) = —,——___& 1 S(tog t) ; 
e 2 eC et Gireye ies) | meen ((ileem ae i = IL) 
1¢F 
E(t) = d(log t) +e | (o(log u))* du 
O S aT olor u))+?((log oy il) 
“al 
G on ee, e ; 
(t) =e 7 t - $ (1 - e((log t) - 2)) 
ok 
+(e *t +1) (1 - $((log t) - 1)) 
sh i 
-Fet(e*t+ 2)(1 - (log t)) +$- 1, 
ope 
H(t) = (e “t + fe) @(log t) - @((log t) - 1). 


for this case the exact form of M(t) cannot be found. 





III. NUMERICAL RESULTS 


Numerical calculations are presented for each bound and 
approximation and all four distributions in Section jit.) ke- 
sults are tabulated below and are graphed and comparisons 
are made with M(t). For the log-normal case M(t) was simu- 


lated for comparison with approximations. 
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FIGURE 9 : M(t) and H(t) for the hyperexponential distribution. 
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